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ABSTRACT 

The  winter-time  transport  of  geostrophically  computed 
angular  momentum  at  500  mb  for  the  period  1  December  1955 
to  28  February  1950  is  obtained  for  Sault  Ste,  Marie, 
Michigan.  This  transport  is  subjected  to  a  cospectral 
analysis  to  compute  the  contribution  of  harmonic  waves  with 
periods  1-50  days  to  the  total  angular  momentum  transport. 
Computations  are  made  for  tv;o  other  stations,  one  east,  and 
one  south  of  Sault  Ste.  Marie  for  comparison  of  transport  as 
related  to  general  circulation  processes. 

The  author  is  de^eply  indebted  to  Professor  F.  L.  Martin, 
Department  of  Meteorology,  U.S.  Naval  Postgraduate  School, 
for  his  encouragement  and  invaluable  aid  in  the  preparation 
of  this  manuscript. 


TABLE  OF  CONTENTS 


Section  Title 

1.  Introduction 

2.  Mathematical  Formulations 

5.  General  Procedures 

4.  Limitations  and  Sources  of  Error 

5.  Interpretation  of  Cospectra 

6.  Conclusions 

7.  Bibliography 

8 .  Appendix  I 


Page 

1 

2 

8 

10 

15 

19 

23 

24 


111 


LIST  OF  ILLUSTRATIONS 


Figure 

Page 

1. 

Illustration  of  X/f  test  of  significance  at  the 
95%  level  of  belief 

12 

2. 

Illustration  of  the  effects  of  sampling  vari¬ 
ations  in  giving  an  aliased  spectrum 

12 

3. 

The  cospectrum  of  the  mean  relative  angular 
momentum  transport  at  the  500-mb  level  over 

Sault  Ste.  Marie,  Michigan,  for  the  period 

1  December  1955  through  28  Februar;y  1956 

20 

4. 

The  cospectrum  of  the  mean  relative  angular 
momentum  transport  at  the  500-mb  level  over 
Caribou,  Maine,  for  the  period  1  December  1955 
through  28  February  1956 

21 

5. 

The  cospectrum  of  the  mean  relative  angular 
momentum  transport  at  the  500-mb  level  over 
Caribou,  Maine,  for  the  period  1  December  195^ 
through  28  February  1955 

21 

6. 

The  cospectrum  of  the  mean  relative  angulax 
momentum  transport  at  the  500-mb  level  over 
Athens,  Georgia, for  the  period  1  December  1955 
through  28  February  1956 

22 

iv 


TABLE  OP  SYMBOLS 


Symbol  Definition 


A\ 

absolute  angular  momentum  per  \init  mass 

U 

zonal  component  of  the  wind 

AT 

meridional  component  of  the  wind 

linear  distance  along  a  parallel  of  latitude 

linear  distance  along  a  meridian 

3 

acceleration, due  to  gravity 

r 

radius  of  the  earth 

t 

density  of  air 

f 

Coriolis  parameter 

grid  length  along  a  parallel  of  latitude 

AY 

grid  length  along  a  meridian 

N 

total  observations 

^3 

geostrophic  wind  component  along  a  parallel  of 
latitude 

geostrophic  wind  component  along  a  meridian 

increment  of  contour  height  along  a  parallel  of 
latitude 

A 

increment  of  contour  height  along  a  meridian 

increment  of  pressure  equal  to  one  millibar 

V 


Symbol  Definition 


A  4 

increment  equal  to  one  centimeter  of  length  along  a 
parallel  of  latitude 

Sir) 

mean  cross  covariance  of  two  continuous  functions 
of  time 

T 

lag  applied  to  a  continuous  function 

UH) 

time  series  of  the  wind  component  along  a  parallel 
of  latitude 

time  series  of  the  wind  component  along  a  meridian 

r 

averaging  period 

Fourier  transform  of  U(t) 

3>(^) 

complex  conjugate  of  the  Fourier  transform  of  Ar(t) 

n 

harmonic  niamber 

c„ 

cospectral  estimate 

T 

fundamental  period 

CO 

angular  velocity 

m 

maximum  lag 

At 

increment  of  time 

A 

integer  of  lag 

mean  cross  covariance  of  two  discrete  functions  of 
time 

Vi 


Symbol 


Definition 


frequency 

chi-square  value  at  the  5%  level 

P  mean  cospectral  value  between  peak  and  trough  of 
the  cospectrum 

V  degrees  of  freedom 

ffn  Nyquist  frequency 


vii 


1.  Introduction 

Widger  [lO]  in  his  study  of  the  flow  of  absolute 

angular  momentim  in  the  atmosphere  stated: 

Although  even  a  complete  knowledge  of  the  angular 
momentum  transfer  in  the  atmosphere  cannot  by  itself 
furnish  a  solution  of  the  problem  of  the  general 
circulation,  any  proposed  scheme  for  this  circulation 
should  include  a  means  for  securing  the  angular 
momentum  flows  which  are  observed.  For  this  reason, 
it  is  essential  that  all  possible  observational  knowl¬ 
edge  concerning' these  processes  be  obtained. 

Many  authors  in  recent  years,  Widger ,• Starr ,  and 

Mintz  [lO,  7,  5]  for  example,  have  provided  studies  in 

connection  with  the  relationship  of  the  absolute  angular 

momentum  transport  to  the  general  circulation,  or  the  global 

balance  of  angular  momentiim.  In  general  these  studies 

provide  quantitative  findings  of  the  spatial  distribution  of 

the  total  momentum  transport  and  establish  the  fact  that 

there  is  a  mean  poleward  transport. 

This  paper,  in  hopes  of  contributing  to  the  knowledge 

of  angular  momentum  transport,  has  as  its  purpose  a  cospect- 

tral  analysis  of  the  angular  momentum  transport  for  waves  of 

periods  1-30  days.  Due  to  the  limited  time  available  for 

the  study,  it  was  confined  to  the  500-mb  level  at  a  station 

in  the  mid-latitudes.  The  primary  station  employed  was 

Sault  Ste.  Marie,  Michigan. 

Several  other  cospectra  are  included  for  other 

stations  for  the  purpose  of  comparison.  However,  the 

spectrum  for  these  comparative  stations  are  for  periods 

2-30  days. 
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2.  Mathematical  Formulations 

Widger  [lo]  has  shown  that  the  rate  at  which  the  ab¬ 
solute  angular  momentum  of  the  zonal  wind  is  increasing  in 
a  portion  of  the  atmosphere  bounded  by  the  earth's  surface, 
the  top  of  the  atmosphere  and  extending  from  the  pole  to  a 
given  latitude  may  be  expressed  as 

J  +  Jprdd-  f  dS 

where  dt/is  any  element  of  volume  within  the  polar  cap  just 
described, d 5  an  element  of  the  surface  of  the  volume,  c/dT 
the  projection  of  an  element  of  the  earth's  surface  area  on 
the  meridional  plane  passing  through  a  given  point,  p  is 
the  pressure,  r  the  radius  of  the  earth,  and  the  east¬ 
ward  frictional  stress  acts  only  on  the  air  at  the  earth's 
surface . 

In  the  first  term  on  the  right  \J^  Is  the  meridional 
component  /v*  of  the  velocity  into  the  polar  cap.  Since  the 
vertical  velocities  at  the  earth' s  surface  and  top  of  the 
atmosphere  may  both  be  considered  zero,  this  quantity  as  a 
whole  represents  the  rate  at  which  angular  momentum  is 
brought  into  the  region  by  air  motion  across  the  boundary. 
The  second  and  third  terms  on  the  right  represent,  respect¬ 
ively,  the  rate  at  which  angular  momentum  is  changed  due  to 
the  torque  exerted  by  pressure  differentials  across  mountain 
ranges  and  the  westv/ard  torque , jrT^,  ,  exerted  by  frictional 
drag  at  the  earth's  surface.  IIov>rever,  since  this  study 
will  consider  only  the  transport  through  a  unit  area  of  the 
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latitude  wall  at  the  500“nit)  surface,  the  terms  involving 
pressure  and  friction  torques  may  he  neglected.  No  attempt 
is  made  herein  to  determine  the  left  side  of  equation  (1) 
from  the  computed  transport  at  ’^00  mb. 

The  absolute  angular  momentiim  of  a  gram  of  air  at 
latitude  ^  may  be  written 


/^1  ~  M  r  C^S  -/i  *  >  r, / 


(2) 


Furthermore ,  since  V„-Ar  and  the  element  of  surface  may  be 


taken  as  an  element  of  the  latitude  wall ,  JS Jz  >  talcing 
the  first  term  on  the,  right  in  equation  (1),  the  transport 
of  absolute  angular  momentum  may  be  v/ritten  as  the  sum  of 
tv/o  terms. 


(5) 


equal  to  the  rate  of  advection  of  "relative"  angular 
momentum,  and  the  iX* momentum  defined  as 


(^) 


v/hich  is  equal  to  the  rate  of  advection  of  angular  momentum 
due  to  the  earth's  rotation.  Since  /L'momentim  does  not 
lend  itself  to  a  cospectral  analysis  and  will  not  be  com¬ 
puted  in  this  study,  the  chief  emphasis  is  placed  upon  the 
average  angular  momentum  transport  through  a  unit  area  and 
per  unit  time  given  by  equation  (3).  Therefore,  after 
applying  the  geostrophic  approximation  for  U  and  AT  in 
equation  (3)  and  dividing  by  N,  one  obtains  the  average 
angular  momentum  transport 
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Ap 


(5) 


v/here  the  hydrostatic  equation  has  also  been  used  in  (5)* 

At  first  glance  equation  (5)  appears  to  he  a  summation 


of  N  terms  around  the  latitude  circle  at  a  distance  A  X 


apart,  and  thus  a  representation  of  equation  (5).  However, 
a  slightly  different  interpretation  was  placed  upon  equa- 
tipn  (5)»  in  that  AX  represents  a  grid  distance  of  500 
nautical  miles  centered  on  the  station  for  v/hich  the 
computation  is  being  made,  and  A  2^  and  are  the  contour- 
height  increments  defining  the  mean  geostrophic  “wind  com¬ 
ponents 


^  -  JL  4,g^ 
9  '  f  AX 


(6) 


and  /vz  :: 


in  the  neighborhood  of  the  station.  In  addition, Ay  was 
taken  equal • to  AX  in  equation  (6).  Thus  the  interpretation 
of  equation  (5)  actually  employed  here  was  to  regard  it  as 
a  time-average,  over  N  equal  time  intervals,  of  values  of 
relative  angular  momentum  transport  centered  at  the  station. 
The  vertical  increment  for  convenience  was  taken  as  Aj?  equal 
to  one  millibar  and  on  the  right  side  of  equation  (5) 
was  scaled  down  to  one  centimeter  instead  of  the  equivalent 
length  of  300  nautical  miles  used  for  the  grid.  The  units 
of  angular  momentum  transport  in  this  study  will  then  be 

Equation  (5)  shov;s  that  the  geostrophic  relative  mom¬ 
entum  transport  depends  upon  the  cross  covariance  betv;een 
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Uj  andA^^.  It  therefore  represents  the  average  geostrophic 
relative  angular  momentum  transported  across  a  section  of 
the  latitude  v;all  with  grid  length  ^/iz  /cm  and  vertical 
thickness  Af^'=  /  m  6  • 

The  mean  cross  covariance  S(t)  of  two  stationary, 
Gaussian  time-series  U(t)  and  at  lagT,  nay  he  ex¬ 

pressed  hy  the  convolution  integral 

^(r)-  'p  M(t) /or  t)  d  t  (7) 

0 

Frequently  in  meteorological  applications,  S(r)  is 
"evenized"  by  formings  the  quantity 


5rr)-f  Sc-t)  (8) 

2 

which  is  independent  of  the  sign  of  lag T .  Henceforth,  it 
will  be  understood  that  the  mean  cross  covariance  has  been 
evenized. 

The  cross-spectrum  (or  simply  cospectrum)  Ch  of  the 
tv/o  series  C((t)  and  is  defined  as  that  part  of  the  mean 
cross  covariance  associated  v;ith  a  unit  v/ave  band  centered 
at  wave  number  h  in  both  series.  In  symbols  may  then 
be  defined  as 


c 


L  I  A' 

p  o? 


(9) 


v/here  P  is  the  averaging  period  and  9^(A)  and  are 

the  Fourier  transforms  of  u(t)  and  respectively,  de¬ 

fined  as  in 


W- 


L  I  m 
T 


-  -X  w-t 

e 


d  t 


«  2Tr^ 
Oi-  —— 


(10) 
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The  circiimflex  indicates  the  complex  conjugate. 

If  one  now  employs  the  definition  of  cross  covari¬ 
ance  5cr)  for  an  evenized  function,  as  contained  in  equation 
(7)»  the  Faltung  theorem  leads  to  the  Fourier  integral 


expression 

5(r)-L /Ywcos  (11) 

SO  that  -SCt')  and  Ccw)  are  spoken  of  as  Fourier-transform 
pairs.  Equation  (11)  may  be  inverted  to  solve  for  Ch  in 
terms  of  Also,  the  resulting  inversion  can  be  spec¬ 

ialized  to  the  case  of  a  finite,  but  long,  averaging  period 
r  t  giving 


rr 


^  J  S(r)  cos 


I  ry>a  A 


(12) 


A  proof  of  equation  (12)  is  given,  from  first  principles, 
in  Appendix  I . 

V/hen  discrete  data  are  employed  at  time  intervals 
of  A.'t  ,  using  a  maximum  lag  of  the  maximum  period, 
or  fundamental,  that  can  be  resolved  is 


-  2.  hn  (13) 

In  equation  (12),  since  T  appears  as  a  summand,  T  will  be 
written  as 


r  =  At 


^  1, . t  (14) 


with  -k.  taking  on  each  of  the  discrete  values  shown  in 
(14).  Since  the  averaging  period  of  equations  (12)  is 
equation  (12)  leads  to 
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(15) 


^  y\ 


[ 


rvc-'V’s^]^z  ^ 


% 


CaS 


17  h-fe  ^ 
J 


/  n  ^  n) 

In  the  case  of  C<j  and  ,  the  resulting  coefficients  are 

divided  by  tv/o,  since  only  one-half  bandwidth  is  available 
at  n~0  and  in~m  ,  Combining  equations  (5)»  (7)i  and  (8), 
the  cross  covariance  for  A  lags  is 


C  r-3r.c.as..^ _ 

'-^4,  2f"CAxf(N->0 


A2,  .  o  AZ  , 


A.flziP(16) 


’  KV) 

sinoe  the  averaging  period  of  equation  (7)  is  now  P-(A/- 
The  factor  two  in  the  denominator  of  equation  (16)  results 
from  evenizing. 

The  cospectral  estimates  offered  by  equation  (15) 
are  termed  "raw  estimates".  Blackman  and  Tucey  Cl]  have 
suggested  several  smoothing  procedures,  the  simplest  of 
which  involves  application  of  a  moving-weighted  average  def-’ 
ined  by 


Cih  -  .'2>5*Cn-i  '^.S  + 


2SC„ 


4-1 


(17) 


In  the  case  of  there  is  no  higher  value  to  include  in 

the  average  and  the  \Treighted  average  is  modified  so  that 

Cm  ^  ^  CjM-i  Cm  (IS) 
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3.  General  Procedures 

The  initial  step  required  establishing  what  data 
were  to  be  used  in  making  computations.  Recalling  that  the 
computed  periods  are  determined  by  the  relation  , 

where  ho  is  the  maximiim  ninnber  of  lags  at  discrete  time 
intervals  At  ,  it  is  clear  that  a  half-day  lag  interval  re¬ 
quires  50  lags  to  obtain  a  spectrum  for  periods  1  -  30  days. 
On  the  other  hand,  a  one-day  sample  and  lag  interval  re¬ 
quires  only  15  lags  to  obtain  a  spectr\im  for  periods  2-30 
days.  Hov;ever,  Blackman  and  Tukey^l^  state  that  the  best 
results  are  obtained  when  the  longest  lag  is  about  5  ~  10 
per  cent  of  the  length  of  the  record.  This  would  require 
that  the  record,  if  possible,  consist  of  300  -  600  observ¬ 
ations  when  using  the  500-mb  charts  produced  at  O3OO  GMT 
and  1500  GMT  daily,  or  half  this  number  when  using  only  one 
series  of  charts.  In  order  to  obtain  this  number  of 
observations,  it  v;ould  require  a  minimum  of  5  months  data 
and  the  spectrum  would  be  influenced  by  seasonal  regimes. 

In  order  to  reduce  the  data-processing  load  and  inter-annual 
influences  as  much  as  possible,  it  was  decided  to  use  the 
observations  from  the  three  v/inter  months  of  December  1955 1 
January  1956,  and  Pebruaiy  1956.  Thus  the  longest  lag  was 
only  17  per  cent  of  the  length  of  the  record  for  this 
particular  season.  However,  no  claim  is  made  that  the 
cospectra  obtained  are  necessarily  representative  of  those 
for  other  winter  seasons. 

Data  were  taken  from  two  sources,  depending  on  the 
station  being  considered.  Since  heights  of  both  the  O3OO 
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GIIT  and  15OO  GMT  500-mb  level  v;ere  required  at  Sault  Ste. 

Marie,  Michigan,  the  primary  station  studied  in  this  paper ^ 
facsimile  charts  transmitted  by  the  U.  S.  Weather  Bureau 
v;ere  utilized.  In  the  case  of  the  other  stations  used  for 
comparison,  the  smallest  period  computed  was  only  two  days 
and  the  data  used  were  the  15OO  GMT  charts  published  in 
Daily  Series  Synoptic  Weather  Maps  [8]. 

From  the  above-mentioned  charts,  the  geopotential 
height  differences  were  obtained  with  a  5*0  degree  (lati¬ 
tude)  grid.  Assuming  that  normal  synoptic  smoothing  \iras 
applied  in  preparation  of  the  charts,  a  grid  of  this  dimen¬ 
sion  should  not  filter  out  significant  short  waves  that 
have  not  already  been  removed  by  synoptic  smoothing. 

In  five  cases,  where  the  O5OO  GMT  facsimile  charts 
were  missing,  geopotential  height  differences  were  computed 
from  the  actual  observed  wind  listed  in  Northern  Hemisphere 
Data  Tabulations  L9j»  The  observed  winds  were  reduced  to 
geostrophic  values  by  subjectively  compensating  for  the 
degree  of  cyclonic  or  anticyclonic  curvature  associated  with 
the  flow. 

The  smoothed  cospectra  of  geostrophic  relative 
angular  momentiam  transport  were  obtained  by  using  these 
data  in  equations  (15) >  (16) »  (17),  and  (18),  and  are  graph¬ 
ically  displayed  in  figs.  5  through  6.  All  computations 
except  for  the  final  smoothing,  utilizing  equations  (17),  (18), 
were  accomplished  on  the  GDC  -  1604  digital  computer. 
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4.  Limitations  and  Sources  of  Error 

The  computational  results  of  this  study  have  been  based 
on  discrete  data  using  equal  intervals  Atr  .  In  one  series 
of  computations,  indicated  by  the  solid  line  in  fig.  3?  the 
data  interval  was  reduced  to  12  hours  in  order  to  determine 
the  effects  of  sampling  fluctuations  and  of  atmospheric 
oscillations.  In  all  other  cospectral  analyses,  the  data 
interval  was  24  hours.  Throughout  the  study  the  maximum  lag 
was  15  days. 

The  period  associated  with  the  smoothed  cospectral 
estimate  is  with  h  =  2^  *  *  ’  •  •  ^  ho  .  On  the  other  hand 

the  central  frequency  of  the  frequency  band  is 


~  2m At  . 

and  the  highest  frequency  resolvable  occurs  for  h = m  ;  this 
frequency  is  called  the  Nyquist  frequency  fm  .  The 
cospectral  estimate  corresponds  to  the  frequency  band 


(h-  54^2 m At  to  (h-f  j4)y2>  mA  t 

It  should  be  noted  that  it  is  difficult  to  resolve 
frequencies  contributing  to  the  co spectrum  at  the  longer 
periods.  For  example  when  h  -  /  ,  the  frequency  band  spans 
the  range  Vao  cycle  /  day  to  cycle  /  day. 

In  order  to  assess  the  reality  of  sharp  peaks  or 
troughs,  Panofsky's  criterion  for  X/y  distribution  of  the 
power  spectrum  has  been  followed.  Strictly  speaking  a  much 
more  rigorous  method  of  assessing  significance  would  have 
been  achieved  by  testing  the  cospectral  correlation 
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function  defined  by 


Here  C>^  is  the  cospectral  estimate  for  the  harmonic, 

and  Ph,^  and  P ^n^^are  the  power  spectral  estimates  for  UL 


and  A/~  at  the  same  harmonic.  Hov;ever,  the  power  spectra 
computations  were  not  available  in  time  for  the  completion  of 
this  aspect  of  the  work.  Hence,  the  more  qualitative  ap¬ 
proach  indicated  in  fig.  1  was  attempted  with  some  misgiv¬ 
ings.  It  must  be  emphasized  that  v/hen  the  term  "95% 
confidence  limit"  is  used,  it  is  only  by  analogy  with  corres¬ 
ponding  phraseology  to  be  found  in  Panofsky  and  Brier  [^4  , 
p.  145]  .  The  actual  confidence  may  be  at  a  somev/hat  lo^^^er 


level. 


distribution  to  the  cospectrum,  a 


Applying  the 


pealc  or  trough  could  not  have  happened  randomly  (at  the  95% 
level  of  belief)  provided  the  peak  or  trough  in  the  cospectral 


estimate  lies  outside  the  range  defined  as 


(19) 


where  r  is  the  mean  cospectral  intensity  between  the  peak 
and  trough,  X"(^s/v  chi-square  value  at  the  .05  level 

for  V  degrees  of  freedom  V  =  2  “  V-LhlZ)  .  Accord¬ 

ing  to  Panofsky  and  Brier  |^4,  p.  
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PERIOD 


Pig.  1.  Illustration  of  the  test  of 

significance  at  the  95%  level  of  belief. 


A  second  source  of  error  which  may  considerably  bias 
the  accuracy  of  cospectral  estimates  is  a  phenomenon  called 
"aliasing".  This  phenomenon  removes  energy  from  higher 
frequencies  than  those  sampled,  and  may  convert  it  to  lower 
frequencies.  Pig.  2  illustrates  how  aliasing  may  occur. 


Pig.  2.  Illustration  of  the  effects  of  sampling 

variations  in  giving  an  aliased  spectrum. 
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Suppose  a  periodic  function,  such  as  that  illustrated 
in  fig.  2,  exists  for  both  uCt)  and  nr(t)  ,  and  consider  the 
period  of  the  oscillation  shown  to  he  2 At  ,  If  lags  are 
talcen  at  intervals  of  3  At  ,  one  is  investigating  the 
cospectral  estimate  at  a  period  of  GAt  .  If  the  sampling 
times  corresponded  to  points  A,  B,  and  C,  one  obtains  a 
rather  large  "aliased"  cospectral  estimate  for  the  period  G At. 
Hov/ever,  if  the  sampling  points  had  been  E,  F,  and  G,  the 
aliased  cospectral  contribution  at  period  GAt  would  be 
zero.  Note  that  the  oscillation  could  be  a  real  atmos¬ 
pheric  oscillation  or, meteorological  "noise"  of  an  obs¬ 
ervational  character.  The  primary  atmospheric  oscillations 
of  periods  shorter  than  one  day  are  the  semi-diurnal  (12-hr) 
wave  and  a  10.5-hr  wave,  originally  described  by  Pekeris  [5]  • 
Hov;ever,  the  amplitude  of  these  waves  is  quite  small. 

Blackman  and  Tukey  [i]  have  shown  that  the  frequencies 
which  are  the  principal  aliases  of  frequency  f„  are 

,  etc.,  where  ^/zAtls  the  Nyquist  freq¬ 
uency.  Hov^ever,  the  effect  of  aliasing  is  most  prono\inced 
at  frequencies  near  the  Nyquist  frequency.  Griffith  et  al 
[^23  suggest  that'  one  way  of  eliminating  aliasing  from  the 
cospectrum  is  to  delete  the  cospectral  estimates  for  the 
last  20  to  ^0  per  cent  of  the  higher  frequencies.  Hov/ever, 
an  experiment  was  performed,  using  data  for  Sault  Ste.  Karie , 
to  assess  the  extent  of  aliasing. 

Fig.  5  shows  tliree  cospectra  for  Sault  Ste.  Marie 
based  on; 

(a)  the  1500  GMT  reports  at  24-hr  intervals  (curve  A), 
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(b)  the  0300  GMT  reports  at  24-hr  intervals  (curve  B) , 
and 

(c)  both  the  0500  GMT  and  I5OO  GMT  reports  at  12-hr 
intervals  (curve  C). 

It  may  be  noted  from  fig*  3  that  all  cospectral  est¬ 
imates  at  periods  greater  than  2, 73  days  are  in  fairly  close 
agreement.  However,  curves  A  and  B  deviate  rather  markedly 
at  their  Kyquist  frequency  oC  Vz  cycle  /  day. 

In  the  case  of  curve  C  there  should  be  little  or  no 
aliasing  at  periods  greater  than  tv;o  days  according  to 
Griffith's  criterion,  since  50  per  cent  of  the  spectrum 
lies  between  the  periods  one  and  two  days.  On  the  other 
hand,  at  all  other  stations  it  vdll  be  assumed,  henceforth, 
that  all  cospectral  estimates  at  periods  2.73  days  or  larger 
are  relatively  unaffected  by  aliasing. 
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5.  Interpretation  of  Cospectra 

Figs.  3  through.  6  present  the  six  co spectra  computed 
for  this  study.  In  these  figures  the  ordinate  is  the  co¬ 
spectrum  and  the  abscissa,  plotted  for  convenience  on  a 
logarithmic  scale,  is  the  period. 

In  assessing  the  reality  of  the  various  peaks  observed 
in  fig.  3»  the  procedure  embodied  in  equation  (19)  vrill  be 
applied  only  to  curve  C,  as  the  other  two  curves  do  not 
deviate  significantly  except  at  a  period  of  3*75  days,  which 
will  be  discussed  separately.  In  order  to  compute  a  mean 
cospectral  estimate  f®r  the  peak  at  a  period  of  10  days,  it 
was  assumed  that  the  cosnectral  estimate  at  30  days  is  the 
apex  of  a  trough.  Inspection  of  the  95%  and  5%  fiducial 
limits  established  from  the  mean  cospectral  values  on  each 
side  of  the  peak  and  trough,  indicate  that  the  peaks  and 
trough  associated  with  periods  10  and  6  days,  respectively, 
are  significant.  The  pealcs  and  troughs  have  only  a  5% 
probability  of  randomly  exceeding  these  limits;  therefore 
both  are  considered  to  be  real. 

W 

In  the  case  of  c\irve  A,  the  peak  appears  not  to  be 
significantly  different  when  compared  to  the  fiducial  limits 
computed  for  curve  C.  However,  limits  computed  from  mean 
values  of  curve  A,  although  not  shown  here  in  order  to  avoid 
confusion  with  those  of  curve  C,  show  the  value  at  the  peak 
on  curve  B  is  considerably  greater  than  the  limit  to  the 
left  side  and  slightly  greater  than  the  limit  to  the  right 
side.  Although  the  pealcs  at  3*75  days  are  not  markedly  sig¬ 
nificant  from  the  limit  on 'the  right  side,  there  is 
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justification  for  stating  these  peaks  seem  to  be  real  when 
one  considers  both  the  left  and  right  side  limits.  Using 
the  significance  test  discussed  above,  the  reality  of  any 
minor  peaks  in  fig.  3  at  periods  of  3  days  or  less  is 
doubtful. 

The  curve  in  fig.  4  is  the  average  cospectrum  for 
Caribou  during  the  period  1  December  1955  bo  28  February 
1956.  A  one-day  interval  between  data  points  was  used  to 
compute  this  cospectrum.  An  inspection  of  the  fiducial 
limits  shown  in  fig.  4  indicates  that  the  peak  at  6  days  is 
significantly  different  from  the  limits  on  both  the  left  and 
right  side  of  the  peak,  and  it  is  therefore  considered  to  be 
real.  As  to  the  small  peak  at  approximately  5*3  days,  the 
confidence  interval  shown  to  the  left  is  quite  representative 
for  this  peak.  It  can  therefore  be  concluded  that  this  peak 
(at  3*5  days)  is  probably  not  real,  as  it  does  not  lie  beyond 
the  confidence  range.  However,  it  can  not  be  inferred  from 
this  that  the  peak  is  a  result  of  sampling  variations,  since 
the  test  does  not  with  certainty  disprove  the  reality  of  the 
pealc. 

The  curve  in  fig.  5  is  also  for  Caribou,  but  xiras  com¬ 
puted  from  data  taken  during  the  winter  previous  to  that 
used  in  fig.  4.  The  limits  in  this  case  show  the  trough  and 
peak  at  7*5  and  5  days,  respectively,  to  be  real.  The  trough 
at  3.75  days  is  questionable,  since  it  is  not  significantly 
different  from  the  limit  to  the  right. 

One  of  the  more  striking  features  in  each  of  the 
cospectra  is  the  large  transport  associated  with  the  longest 
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period  and  the  fact  this  transport  has  the  same  sign  as  the 
total  transport. 

In  general  it  is  difficult  to  associate  a  particular 
period  of  the  cospectrum  with  a  specific  space  wave,  since 
the  period  is  a  function  of  both  the  wavelength  and  speed 
of  the  space  wave.  It,  thereby,  seems  possible  for  a  trans¬ 
port  at  a  specific  period  to  have  contributions  from  both 
long  and  short  space  v^aves.  This  would  be  the  case  if  the 
ratio  of  their  wavelengths  is  equal  to  the  ratio  of  their 
speeds.  However,  assuming  that  the  Rossby  wave  equation 
gives  a  valid  relationship  between  the  wave  speed  and  wave¬ 
length  for  any  given  mean  zonal  wind,  it  may  be  expected 
that  long  space  v/aves  move  eastward  slowly  and  short  ones 
more  rapidly.  In  view  of  this  assumption,  the  transport 
associated  with  a  period  of  thirty  days  may  be  attributed 
to  the  contributions  from  long  space  waves. 

Purthermore ,  large  values  of  transport  by  the  longest 
period  were  to  be  expected,  for  an  examination  of  the 
monthly  mean  500-mb  charts  of  December  1955  through  1956 
showed  Sault  Ste.  Marie  and  Caribou  were  on  the  V7e stern 
side  of  a  mean  long  wave  trough,  hence  the  net  transport 
would  most  likely  be  southward  (negative).  However,  for  the 
second  cospectrum  at  Caribou  (fig.  5)>  the  period  December 
195^  through  February  1955  was  selected  by  ascertaining 
through  monthly  mean  500-mb  charts  that  Caribou  was  east  of 
the  trough  and  a  northward  (positive)  transport  could  be 
expected.  In  the  case  of  Athens,  Georgia,  the  main  reason 
for  selecting  this  station  to  provide  a  comparative 
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cospectrum  is  the  fact  it  is  almost  due  south  of  Sault  Ste. 
riarie.  One  could  expect  these  two  stations  to  have  a  net 
transport  of  the  same  sign,  but  examination  of  the  monthly 
mean  500-mb  charts  show  the  mean  long  wave  trough  to  have  a 
southv/est  to  northeast  tilt  of  sufficient  magnitude  to  place 
the  trough  west  of  Athens.  As  a  result  the  net  transport  for 
Athens  v/as  northward,  while  that  of  Sault  Ste.  Marie  was 
southward.  The  indicated  convergence  of  momentum  transports 
should  then  act  in  such  a  v;ay  as  to  increase  the  mean  zonal 
v;esterlies.  Note  also  that  the  only  pronounced  peak  in  the 
cospectrum  at  Athens  occurs  at  a  period  of  50  days,  and  that 
the  cospectral  values  diminish  v;ith  decreasing  period.  This 
is  probably  due  to  the  fact  that  short  period  waves  do  not 
penetrate  appreciably  into  this  latitude  belt. 

Interpretation  of  significant  transports  associated 
with  shorter  periods  requires  reference  to  the  earlier 
statement,  that  short  v;aves  are  likely  to  move  rapidly.  In 
viev;  of  this  concept  the  momentum  transport  evident  near 
periods  of  four  days  is  attributed  to  waves  of  approximate 
space-wave  numbers  6  to  10,  which  Saltzmann  and  Fleisher 
[6]  have  called  "cyclone  v^aves".  However,  the  moment\im 
transport  peaks  at  periods  near  10  days  are  presumably  due 
to  long  waves  of  space-vrave  numbers  1  to  5» 
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6.  Conclusions 

It  nay  be  concluded,  tentatively,  from  the  cospectra 
presented  in  figs.  3  through  6  and  the  preceding  discussion 
that : 

(a)  Haximm  transport  of  geostrophic  relative  angular 
noHentun  is  accomplished  by  disturbances  xirith  the  longer 
periods, 

(b)  The  pealcs  with  periods  near  10  days  are  related 
to  long  space  waves. 

(c)  The  direction  of  the  net  transport  is  hi^ly  de¬ 
pendent  on  the  position  of  the  mean  long  wave  trough 
relative  to  the  observing  station. 

(d)  Cyclone  v;aves  of  space-wave  numbers  6  to  10  may 
contribute  to  the  maximum  momentum  transports  associated 
x;ith  periods  near  4  days. 

(e)  The  maximum  transports  at  periods  of  50  days  may 
be  attributed  to  the  circulation  index  cycle. 


19 


ro 

o 


Fig.  5  The  cospectrum  of  the  mean  relative  angular  momentum  transport  at  the  500-ni^  level 
over  Sault  Ste.  Marie,  Michigan  for  the  period  1  December  1955  through  28  February 
1956.  The  aiixiliary  graph  has  the  ordinate  enlarged  100  times  for  periods  1  to  2 
days. 


Fig.  4  Th,e  cospectrum  of  the  mean  relative  angular  momentum 
transport  at  the  500-mb  level  over  Caribou,  Maine  for 
the  period  1  December  1955  through  28  February  1956. 


Fig. 


5  The  cospectrum  of-  the  mean  relative  angular  momentum 
transport  at  the  500-mb  level  over  Caribou ,  T'Taine 
for  the  period  1  December  195^  through  28  February  1955o 
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Fig.  6  The  cospectrum  of  *the  mean  relative  angular  momentum 
transport  at  the  500-mb  level  over  Athens,  Georgia 
for  the  period  1  December  1955  through  28  February  1956. 
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APPENDIX  I 


Derivation  of  the  Cospectral  Function 

The  cross  covariance  at  lagT,  of  two  functions  of 
time,  is  given  by 


S(x)  if,  U(t)  .  ir('t  #■  “)  d  ■= 


^  S,,)]  [|8„  t 

1 1 /?„ ^ -  => )  t)]  J 

-  rV-^t/Z  ^ 4  .2«^ 

>  »o  h  ^  ^  '  »  -» 

4  cc»5  ej ~~  1  ^ 


■i^Jj 

The  first  cosine  term  integrates  to  zero  for  all  cases 


m,n  ^  0.  The  second  cosine  term  integrates  to  zero  for 

all  cases  m  n.  When  m  =  n,  integration  gives 

r  r  T 

{f.- 

t)  "^O  ^ 

—  *  >« 

5^4,1  (t- 

b-a  ' 

The  evenized  cross  covariance  equals  one-half  the  sum 
of  the  cross  covariance  at  T  and-T  , 

* 

\  -  5(<-r)  -  L  V /^/,  By,  /  , /,  ,,  n  zphX  ^ 

i^Cr)  2^  ^  7  y-  C  O  5  (,  /^  S/i )  c  oS  y  H  .  Zi  h  -  /,  2,  ‘  •  • ' 

h-0 


Defining  the  unsmoothed  cospectrinn  as 


Cos()^  ) 


24 


the  Fourier  transform  of  C(h)  is  obtained  as 


S'w-  r  Ccn)  -^-  dh 

Multiplying  both  sides  by  C(3S~^  J  V  and  integrating  from 
0  to  T  gives 

S{r)ca5  (t  ■*"  dT 

4iY  hT 

For  all  n  /  0,CoS - —  integrates  to  zero,  giving 

C(n)-  I  S(r)  COS—^-  dr 

v/here  P  is  the  integration  limit  and  T  is  the  fundamental 
peri od. 
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